TricLxial projected shell model study of 7-vibrational bands in even-even Er isotopes 



00 

o 
o 

(N 



> 

p 

O 

00 

o 



X 



J. A. Sheikh\ G. H. Bhat^, Y. Sun^^^, G. B. Vakil\ and R. Palif* 
^Department of Physics, University of Kashmir, Srinagar, 190 006, India 
^Department of Physics, Shanghai Jiao Tong University, Shanghai 200240, People's Republic of China 
^ Joint Institute for Nuclear Astrophysics, University of Notre Dame, Notre Dame, Indiana 46556, USA 
'^Tata Institute of Fundamental Research, Colaba, Mumbai, 400 005, India 

We expand the triaxial projected shell model basis to include triaxially-deformed multi- 
quasiparticle states. This allows us to study the yrast and 7-vibrational bands up to high spins 
for both 7-soft and well- deformed nuclei. As the first application, a systematic study of the high- 
spin states in Er-isotopes is performed. The calculated yrast and 7-bands are compared with the 
known experimental data, and it is shown that the agreement between theory and experiment is 
quite satisfactory. The calculation leads to predictions for bands based on one- and two-7 phonon 
where current data are still sparse. It is observed that 7-bands for neutron-deficient isotopes of 
^^^Er and ^^^^Er are close to the yrast band, and further these bands are predicted to be nearly 
degenerate for high-spin states. 

PACS numbers; 21.60.Cs, 21.10.Hw, 21.10.Ky, 27.50.-|-e 



I. INTRODUCTION 

Recent experimental advances in nuclear spectroscopic 
techniques following Coulomb excitations, in-elastic neu- 
tron scattering, and thermal neutron capture have made 
it possible to carry out a detailed investigation of 7- 
vibrational bands in atomic nuclei [J, [Ij 3 ■ These bands 
are observed in both spherical and as well as in deformed 
nuclei. In spherical nuclei, the vibrational modes are well 
described using the harmonic phonon model [1, Q . Al- 
though exact harmonic motion has never been observed, 
there are numerous examples of nuclei exhibiting near 
harmonic vibrational motion. As a matter of fact, one- 
and two-phonon excitations have been reported in a large 
class of spherical nuclei. In deformed nuclei, vibrational 
motion is possible around the equilibrium of deformed 
shape configuration. The deformed intrinsic shape is pa- 
rameterized in terms of P and 7 deformation variables. 
These parameters are related to the axial and non-axial 
shapes of a deformed nucleus. The one-phonon vibra- 
tional mode in deformed nuclei with no component of 
angular momentum along the symmetry axis {K = 0) is 
called /3-vibration and the vibrational mode with com- 
ponent of angular momentum along the symmetry axis 
{K = 2) is referred to as 7-vibration. The rotational 
bands based on the 7-vibrational state are known as 7- 
bands [1, 0) Hi • One-phonon 7-bands have been observed 
in numerous deformed nuclei in most of the regions of the 
periodic table. There has also been reports on observa- 
tion of two-phonon 7-bands [1, [l3| . 

Several theoretical models have been proposed to study 
7- bands with varying degree of success. The quasipar- 
ticle phonon nuclear model (QPNM) [HI, [l^, which re- 
stricts the basis to, at the most, two phonon states, has 
led to the conclusion that two-phonon collective vibra- 
tional excitations cannot exist in deformed nuclei due 
to the Pauli blocking of important quasiparticle com- 
ponents. On the other hand, the multi-phonon method 
(MPM) embodies an entirely different truncation 



scheme. It employs only a few collective phonons and re- 
stricts the basis to all the corresponding multi-phonon 
states up to eight phonons. This approach predicts that, 
for strongly collective vibrations, two phonon K'^ = 4+ 
excitations should appear at an energy of about 2.6 times 
the energy of the one-phonon K'^ — 2+ state [3, [H] . On 
the other hand, the dynamic deformation model (DDM) 
p^ . which is quite different from the models mentioned 
above, constructs collective potential from a set of de- 
formed single-particle basis states accommodating eight 
major oscillator shells. This model predicts a collective 
K'^ = 4+ at almost 2 MeV. 

AU the above mentioned models (QPNM, MPM, and 
DDM) do not have their wave functions as eigen-states of 
angular momentum. Strictly speaking, these methods do 
not calculate the states of angular momentum, but the 
if-states (K is the projection of angular momentum on 
the intrinsic symmetry axis). To apply these models, one 
has to assume that I K. However, since an intrinsic 
A'-state can generally have its components spread over 
the space of angular momenta of / > AT, the reliability of 
these approaches depends critically on actual situation. 
As pointed out by Soloviev [l2|, it is quit desirable to 
recover the good angular-momentum in the wave func- 
tions. 

Some algebraic models including the extended version 
of the interacting boson (sdg-IBM) [iBl and pseudo- 
symplectic models have also been employed to study 
the 7- excitation modes and predict high collectivity for 
the double 7- vibration (20j . 

Recently, the triaxial projected shell model (TPSM) 
has been employed to describe 7- bands [2l|, |2^. This 
model uses shell model diagonalization approach and in 
this sense, it is similar to the conventional shell model 
approach except that the basis states in the TPSM 
are triaxially deformed rather than spherical. In the 
present version of the model, the intrinsic deformed ba- 
sis is constructed from the triaxial Nilsson potential. 
The good angular momentum states are then obtained 
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through exact three-dimensional angular momentum pro- 
jection technique. In the final stage, the configuration 
mixing is performed by diagonalizing the pairing plus 
quadrupole-quadrupole Hamiltonian in the projected ba- 
sis [H Hi]. The advantage of the TPSM is that it 
describes the deformed single-particle states microscopi- 
cally as in QPNM, MPM, and DDM, but its total many- 
body states are exact eigen states of angular momentum 
operator. Correlations beyond the mean-field are intro- 
duced by mixing the projected configurations. 

It is to be noted that an intrinsic triaxial state in the 
TPSM is a rich superposition of different K-st&tes. For 
instance, the triaxial deformed vacuum state is composed 
of iiT = 0, 2, 4, • • • configurations. The projected bands 
from these K = 0,2 and 4 intrinsic states are the dom- 
inant components of the ground-, 7-, and 27-bands, re- 
spectively [22]. 

In the earlier TPSM analysis for even-even nuclei, the 
shell model space was very restrictive, including only 0- 
quasiparticle (qp) state [2li [H, [H, |2i, 113, [l^. This 
strongly limited the application of the TPSM to the low- 
spin and low-excitation region only. It was not possi- 
ble to study high-spin states because multi-qp config- 
urations will usually become important for states with 
J > 10 in the normally deformed rare-earth nuclei. In the 
present work, the qp-space is enlarged to incorporate the 
two-neutron-qp, two-proton-qp and four-qp configuration 
consisting of two protons plus two neutrons. This large 
qp space is adequate to describe the bands up to second 
bandcrossing [2l|. The purpose of the present work is, 
as a first application of the extended model, to perform 
a detailed investigation of the high-spin band structures, 
in particular 7-bands, of Erbium isotopes ranging from 
mass number A = 156 to 170. In a parallel work [2^, the 
TPSM analysis for odd-odd nuclei in a multi-qp space 
has been performed. 

The manuscript is organized in the following manner: 
in the next section, a brief description of the TPSM 
method is presented. The results of the TPSM study 
are presented and discussed in section III. Finally, the 
work is summarized in section IV. 



II. TRIAXIAL PROJECTED SHELL MODEL 
APPROACH 



In the present work, the TPSM qp basis is extended, 
which consists of projected 0-qp vacuum, 2-proton (2p), 
2- neutron (2n), and 4-qp states, i.e.. 



In Eq. (fit, th e three-dimensional angular-momentum op- 
erator is [301 
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with the rotational operator 
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and I $ > represents the triaxial qp vacuum state. The 
qp basis chosen above are adequate to describe the high- 
spin states upto, say / ~ 24, and in the present analysis 
we shall restrict to this spin regime. The triaxially de- 
formed qp states are generated by the Nilsson Hamilto- 
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Here Hq is the spherical single-particle Hamiltonian, 
which contains a proper spin-orbit force |31j. The pa- 
rameters e and e' describe axial quadrupole and triax- 
ial deformations, respectively. It should be noted that 
for the case of axial-symmetry, the qp vacuum state has 
K = 0, whereas in the present case of triaxial deforma- 
tion, the vacuum state | $ > is a superposition of all the 
possible iiT- values. The allowed values of the if -quantum 
number for a given intrinsic state are obtained through 
the following symmetry consideration. For the symmetry 
operator, S = e~'^'^'^' , we have 

PiiK \^>^PLkS^S\<^>> = e-^^-'^^Pi,^ I $ > , (5) 

where 5 | $ > = g^"'' | $ > , and k characterizes the in- 
trinsic states in Eq. ([1]). For the self-conjugate vacuum or 
0-qp state, k = and, therefore, it follows from the above 
equation that only K = even values are permitted for this 
state. For 2-qp states, the possible values for X-quantum 
number are both even and odd depending on the struc- 
ture of the qp state. For the 2-qp state formed from the 
combination of the normal and the time-reversed states, 
K — and, therefore, only K = even values are permit- 
ted. For the combination of the two normal states, k = 1 
and only K = odd states are permitted. 

As in the earlier projected shell model (PSM) calcu- 
lations, we use the pairing plus quadrupole-quadrupole 
Hamiltonian [2^ 

H = Ho-lxYl QlQ^^ - GmP^P ^GqY, PlP^.- (6) 



The interaction strengths are taken as follows: The 
QQ-force strength x is adjusted such that the physi- 
cal quadrupole deformation e is obtained as a result of 
the self-consistent mean- field HFB calculation f23|. The 
monopole pairing strength Gm is of the standard form 

Gm = [21.24 T 13.86(7V - Z)/A] /A, 
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with "— " for neutrons and "+" for protons, which ap- 
proximately reproduces the observed odd-even mass dif- 
ferences in the rare-earth mass region. This choice of 
Gm is appropriate for the single-particle space employed 
in the PSM, where three major shells are used for each 
type of nucleons {N = 4, 5, 6 for neutrons and = 3, 4, 5 
for protons). The quadrupole pairing strength Gq is as- 
sumed to be proportional to Gm , and the proportionality 
constant being fixed as 0.18. These interaction strengths 
are consistent with those used earlier for the same mass 
region [HI, [H, [H . 

III. RESULTS AND DISCUSSIONS 

The triaxial projected shell model calculations have 
been performed for Er-isotopes ranging from A = 156 
to 170. The deformation parameters (e, e') used in the 
present work are same as those employed in Ref. [23 |. 
It has already been mentioned in section II that in the 
present work the mean-field potential is constructed with 
given input deformation values of e and e'. In a more re- 
alistic calculation, these deformation values for a given 
system are obtained through the variational HFB calcu- 
lations. The chosen values of e for the present calculation 
are those from the measured quadrupole deformations of 
the nuclei as is done in the previous projected shell model 
analysis. The e' values used in the present work are real- 
istic, which correctly reproduce, for example, excitations 
of the 7 band relative to the ground-state [22 |. 

A. Band Diagrams 

Band diagrams can bring valuable information regard- 
ing the underlying physics [2^. These band diagrams 
for the studied Er-isotopes are presented in Figs. 1 to 4 
and depict the results of the projected energies for each 
intrinsic configuration. In the diagrams, the projected 
energies are shown for 0, 2n, 2p and 2p+2n quasiparticle 
configurations. The qp energies for these configurations 
are given in the legend of each figure. As already men- 
tioned in the last section that with the triaxial basis, the 
intrinsic states do not have a well-defined i^-quantum 
number. Each triaxial configuration in Eq.([T]) is a com- 
position of several X-values and bands in Figs. 1 to 4 
are obtained by assigning a given i^T- value in the angular- 
momentum projection oprator. To make the discussion 
easy, we denote a if -state of an i-configuration as {K, i), 
with i = 0, 2n, 2p and 4. For example, the K ^ state 
of 0-qp configuration is marked as (0, 0) and K = 1 oi 
2n-qp configuration as (1, 2n). 

In Figs. 1 to 4, the projected bands associated with the 
0-qp configuration are shown for K ^ 0, 2, and 4, namely 
the (0, 0), (2, 0), and (4, 0) bands. In the literature, these 
K = 0,2, and 4 bands are referred to as ground-state, 7-, 
and 27-bands. The ground-state band has k = and is, 
therefore, comprised of only even-K values. We use the 



same names in the following discussion to be consistent 
with the literature, but stress that in our final results ob- 
tained after diagonalisation, K is not a strictly conserved 
quantum number due to configuration mixing. 

It is evident from Fig. 1 that the (2, 0) bands for ^^^Er 
and ^^^Er lie very close to the (0, 0) bands. This means 
that 7- vibration has low excitation energy in these two 
nuclei. For high-spin states, it is further noted that the 
(0, 0) and (2, 0) band energies become almost degenerate, 
and as a matter fact for 1—16 and above, the energy 
of even-spin states in the (2, 0) band is slightly lower 
than the (0, 0) band. It is a well-known fact that 7- 
bands become lower in energy with increasing triaxility 
and what is also evident from Fig. 1 that they become 
favored with increasing angular-momentum. As can be 
seen from Fig. 1, the (2, 0) bands in ^^^Er and ^^^Er also 
depict pronounced signature splitting with the splitting 
amplitude increasing with spin. The (4, 0) band is close 
to the (2, 0) band for ^'"'^Er and lies at a slightly higher 
excitation energy for ^^^Er. The (4, 0) bands in these 
two isotopes are also noted to have signature splitting 
for higher angular momenta, and the splitting amplitude 
is nearly the same for the (2, 0) and (4, 0) bands. 

In Fig. 1, several representative multi-qp bands, 
namely projected 2- and 4-qp configurations, are also 
plotted. Although the K = 1 2-qp neutron (1, 2n) and 2- 
qp proton (1, 2p) bands are close in energy for low spins, 
but with increasing spin the 2n-qp bands are lower in en- 
ergy than 2p-qp bands due to larger rotational alignment. 
It is noted that neutrons are occupying lii3/2 and pro- 
tons are occupying lhii/2 intruder sub-shells. For each 
of the (1, 2n) and (1, 2p) bands, the projected energies 
are also shown for the corresponding 7-bands with con- 
figurations (3, 2n) and (3, 2p). The (1, 2n) band is noted 
to cross the (2, 0) and the (0, 0) bands at / = 12. It is 
also seen that the (3, 2n) band crosses the (0, 0) band 
at a slightly higher spin value of / = 14. It is interesting 
to note that after the band crossing, the lowest even-spin 
states originate from the (1, 2n) band, whereas the odd- 
spin members are the projected states from the (3, 2n) 
configuration. Finally, the 4-qp (4, 4) configuration lies 
at high excitation energies and does not become yrast, 
at-least, up to the spin values shown in the figure. 

The band diagrams for ^^°Er and ^^^Er are presented 
in Fig. 2. The energy separation between the (0, 0) and (2, 
0) bands is larger as compared to the two lighter isotopes 
in Fig. 1. In the case of ^^^Ei, the (2, 0) band energies 
do come close to the (0, 0) energies for spins / > 12. 
The (1, 2n) band again crosses the (0, 0) band at / = 12 
for ^^"Er and at / = 14 for ^^^Er. The band diagrams 
for ^^^Er and ^^^Er shown in Fig. 3 depict larger energy 
gaps among various bands. The signature splitting of the 
(2, 0) band has considerably reduced. It is further noted 
that 2n-band-crossing is shifted to higher spin values. For 
the case of "'^^'^Er, the band crossing is observed to occur 
at / = 16 and for ^^^Er it occurs at / = 18. The band 
diagrams for ^^^Er and ^™Er shown in Fig. 4 indicate 
that the (2, 0) bands are quite high in excitation energy. 
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The band crossing for these cases is further shifted to 
higher spin values. 



B. Results after Configuration Mixing 

In the second stage of the calculation, the projected 
states obtained above are employed to diagonalize the 
shell model Hamiltonian of Eq. © . It is to be mentioned 
that for the discussion purpose, only the lowest three 
bands from the 0-qp configuration and lowest two bands 
for each other configuration have been shown in band 
diagram. Figs. 1 to 4. However, in the diagonalisation 
of the Hamiltonian, the basis states employed are much 
more, which includes, for example, those K ~ 1,3,5 and 
7 with K = l and = 0, 2, 4, 6 and 8 with k = 0. 

The lowest three bands after the configuration mixing 
are shown in Figs. 5 and 6 and are compared with the 
experimental energies wherever available. Although they 
are of mixed configurations in our model, we still call 
them yrast, 7- and 27-bands to be consistent with the 
literature. It is observed from these two figures that the 
agreement between the calculated and the experimen- 
tal energies for the yrast and 7-bands is quite satisfac- 
tory. For ^YiQ theoretical yrast line depicts two 
slopes and these correspond to the slopes of two crossing 
bands shown in Figs. 1 and 4. This also indicates that the 
interaction between the two crossing bands is small with 
the result that these nuclei shall depict a back-bending 
effect [13]. It is also encouraging to note from Figs. 5 
and 6 that the agreement for the 7-bands is quite good, 
except that for ^^"'Er and ^^''Er, the signature splitting 
at the top of the bands is not reproduced properly. For 
the 27-bands, our calculations agree well with the only 
available data in ^^^Er d and ^^^Er [lO|]. 

There is another notable effect about anharmonicity 
in 7 vibrations. If we regard the 7-bandhead as one 7- 
phonon vibration and the 27-bandhead as two 7-phonon 
vibration, it can be easily seen from Figs. 5 and 6 that the 
vibration is not perfectly harmonic. In fact, in the two 
lightest isotopes, the 7-soft ^^^Er and ^^^Er, the vibration 
is almost harmonic. As the neutron number increases, 
a clear anharmonicity is predicted from our calculation 
and the degree of anharmonicity increases with increasing 
neutron number. 



C. Analysis of Wavefunction 



inant from / = 12 to 16. For / = 18 and onwards, there 
are many configurations with finite values contributing 
to the yrast states. The band diagram of ^^^Er in Fig. 
1 suggests that the 7-band should have the (2, 0) con- 
figuration as the dominant component. This is evident 
from Fig. 7 and it is also noted that (0, 0) is significant 
for the even spin states up to / = 8. The / = 10 state 
is mostly composed of (1, 2n) and for higher spin states 
the (3, 2n) and (3, 2p) configurations are the dominant 
components of the 7 band. The 27 band in Fig. 7 is 
composed of (4, 0) band for the low spin states. 1 — 8 
of this band is predominantly composed of the (1, 2n) 
configuration, but the high spin states are found to have 
quite a complex structure. 

The yrast wavefunction decomposition of ^®^Er, shown 
in the top panel of Fig. 8, indicates that this lowest band 
is predominantly composed of the (0, 0) configuration up 
to / = 12 and there appears to be very small admixtures 
of K — 2 and other configurations. After the bandcross- 
ing at / = 16, the yrast states are dominated by the (1, 
2n) configuration. There is also a significant contribu- 
tion of the (3, 2n) configuration after the bandcrossing. 
The 7-band in Fig. 8 is primarily composed of the (2, 
0) configuration up to / = 11 and above this spin the 
states are a mixture of different configurations. There is 
a clear distinction in the composition of the even- and 
odd-spin states above / = 11. The odd-spin states are 
composed of the (3, 2n) and (2, 0) configurations, and 
the even-spin states are dominated by the (1, 2n) and (0, 
0) structures. The 27-band up to / = 7 is primarily the 
(4, 0) configuration. For 7 = 8 and above, this band is a 
mixture of (1, 2n) and (3, 2n) configurations. 

The wavefunction analysis of ^^"^Er shown in Fig. 9 
indicates that the yrast state, as expected for a well de- 
formed nuclei, is mainly comprised of the (0, 0) config- 
uration. This contribution drops smoothly and, on the 
other hand, the (1, 2n) component increases steadily. For 
/ — 20, it is noted that the (0, 0) and (1, 2n) contribu- 
tions are almost identical and above this spin value, it 
is expected that the (1, 2n) configuration shall dominate 
the yrast states. The 7-band is also noted to have a well 
defined structure of (2, 0) and only for high spin states, 
it is observed that the (1, 2n) and (3, 2n) of the 2n- 
aligned configuration become important. The 27-band is 
dominated mostly by the aligning configurations above 
1 = 7. As is evident from the band diagram of this nu- 
cleus, presented in Fig. 4, that 2n-aligned band is lower 
than the (4, 0) band for most of the spin values. 



In order to probe further the structure of the bands 
presented in Figs. 5 and 6, the wavefunction decompo- 
sition of the yrast, 7- and 27-bands are shown in Figs. 
7, 8 and 9 for ^^^Er, ^^^F,t and "°Er. For other nuclei, 
the wavefunction have smiliar structure and are not pre- 
sented. It is seen from Fig. 7 that the yrast band for 
-'^^^Er is predominantly composed of the (0, 0) configu- 
ration up to / = 10. The (0, 0) contribution suddenly 
drops at / = 10, and (1, 2n) configuration becomes dom- 



IV. SUMMARY AND CONCLUSIONS 

In the present work, the triaxial projected shell model 
approach with extended basis has been employed to 
study the high-spin band structures of the Er-isotopes 
from A = 156 to 170. In this model, the Hamiltonian 
employed consists of pairing plus quadrupole-quadrupole 
interaction. It is known that Nilsson deformed potential 
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is the mean-field of the quadnipole-quadrupole interac- 
tion and this potential is directly used as the Hartree- 
Fock field rather than performing the variational calcu- 
lations. It is, as a matter of fact, quite appropriate to 
use the Nilsson states as a starting basis because the pa- 
rameters of this potential have been fitted to large body 
of experimental data. The parameters of the model are 
the deformation parameters of e and e'. The axial de- 
formation parameter e has been fixed from the observed 
quadrupole deformation of the system as is done in most 
of the projected shell model analysis. The non-axial 
parameter e' was chosen to reproduce the bandhead of 
the 7 band. The pairing strength parameters have been 
determined to reproduce the odd-even mass differences. 
The monopole pairing interaction has been solved in the 
BCS approximation and the qp states generated. In the 
present work, the qp states considered are: 0-qp, 2-qp 
neutron, 2-qp proton, and the 4-qp state of 2-neutron 
plus 2-proton. 

In the second stage of the calculations, the three- 
dimensional angular-momentum projection is performed 
to project out the good angular-momentum states from 
these qp states. These projected states are then used 
as the basis to diagonalise the shell model Hamiltonian 
in the third and the final stage. The salient features of 
results obtained in the present work are: 



1. 7-bands are quite close to the yrast line for the 
neutron-deficient Er-isotopes, in particular, for 
^^^Er and ^"''^Er. It is further evident from the 
present results that these 7-states become further 
lower in energies for high-spin states. As a matter 
of fact, for ^^^Er and '^^^Er, they become lower than 
the ground-state band for I > 14. We propose that 
this is a feature of 7-soft nuclei. 

2. 7-bands are pushed up in energy with increasing 
neutron number, and further the degree of anhar- 
monicity of 7 vibration also increases. 

3. The wavefunction decomposition of the bands 
demonstrates that for neutron deficient Er- 
isotopes, there is a significant mixture of the 7 con- 
figuration in the ground-state band and vice-versa. 
The neutron rich ^^°Er nucleus, on the other hand, 
has the intrinsic structures as expected for a well 
deformed nucleus with the ground-state band com- 
posed of nearly pure K = configuration. 
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FIG. 1: Band diagrams for ise-issgj. jgotopes. The labels (0,0), (2,0), (4,0), (l,2n), (3,2n), (l,2p), (3,2p), (2,4) and (4,4) 
correspond to ground, 7, 27, two neutron-aligned, 7-band on this two neutron-aligned state, two proton-aligned, 7-band on 
two this proton-aligned state, two-neutron plus two-proton aligned band and 7 band built on this four-quasiparticle state. 
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FIG. 3; Band diagrams for Er isotopes. The labels indicate the bands mentioned in the caption of Fig. 1. 
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FIG. 8: Wavefunction decomposition for Er. ax denotes the amplitude of the wavefuction in terms of the projected basis 
states. 
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FIG. 9: Wavefunction decomposition for Er. ax denotes the amplitude of the wavefuction in terms of the projected basis 
states. 



